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O ■ Abstract 

0^ ' We present the first-order coherent exclusive exponentiation (CEEX) scheme, 

r~| ' with the full control over spin polarization for all fermions. In particular it is 

Q-i! applicable to difficult case of narrow resonances. The resulting spin amplitudes 

QhI and the differential distributions are given in a form ready for their implementation 

^ ■ in the Monte Carlo event generator. The initial-final state interferences are under 

control. The way is open to the use of the exact amplitudes for two and more hard 
photons, using Weyl-spinor techniques, without giving up the advantages of the 
^ ■ exclusive exponentiation, of the Yennie-Prautschi-Suura type. 
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1 Introduction: What is the problem? 



The problem addressed in this work is: How to describe, consistently in the process 
e"'"e~ — >■ //, the coherent emission of initial state radiation (ISR) and final state radiation 
(FSR) of soft and hard photons, providing for cancellations of infrared (IR) divergences 
from real and virtual photon emission to infinite perturbative order (exponentiation), at 
the level of completely exclusive multiphoton differential distributions, i.e. in the form 
suitable for implementation in Monte Carlo (MC) event generators? In addition we are 
looking for the solution that is friendliest to narrow s-channel resonances. 

This work is firmly rooted in the work of Yennie, Frautschi and Suura (YFS) on QED 
exponentiation and its further developments in refs. . The present work definitely 
goes beyond the scope of these previous papers - the main difference is the consequent 
use of spin amplitudes in the exponentiation. Our work is close in spirit, although not in 
technical details, to seminal papers of Greco et al. [^,|^ on QED exponentiation for narrow 
resonances. However, it should be stressed that, contrary to refs. 0,0, all our differential 
multiphoton distributions are completely exclusive (important for MC implementation) 
and we do include hard photons completely and systematically. In this context, the 
work of ref. should also be mentioned. It implements QED interferences among e+ 
and e~ fermion lines, the analog of the ISR-FSR interferences, for the first time in the 
exclusive exponentiation. It does not, however, use spin amplitudes for exponentiation 
as consequently as does the present work; it is also rather strongly limited to exact first 
order exponentiation in the YFS framework. It is sort of half-way between the present 
work and the older ones of refs. [0,0. At the technical level, the methods used here for the 
construction of the spin amplitudes are essentially thoseQ of Kleiss and Stirling (KS) 
with the important supplement of ref. |jlO|, providing for total control of complex phases 



and/or fermion spin quantization frames. The MC implementation of the present work 
will soon be available [|lT| and it will replace two MC programs: KORALB |jl2|], where 
fermion spin polarizations are implemented exactly, but there is no exponentiation, and 
KORALZ |jl3[, where exponentiation is included, but the treatment of spin effects is 
simplified^. 

The present work is essential for any present experiments in e'^e~ colliders and future 
e~^e~ and fj^'^fi~ colliders, where the most important new features for data analysis will 
be inclusion of ISR-FSR interferences and (in the next step) the exact matrix element for 
emission of 2 and 3 hard photons, in the presence of many additional soft ones. 



2 Basic KS/GPS spinors and photon polarizations 

The arbitrary massless spinor u\{p) of momentum p and chirality A is defined according to 
KS methods In the following we follow closely the notation of ref. |]10[ (in particular 

^ We have evaluated several techniques based on Weyl-spinor techniques and we concluded that the 
technique of KS is best suited for our needs (exponentiation) . 

^ The ISR-FSR interference is also neglected in KORALZ, in the main mode with the exponentiation 
switched on. 
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we also use C = O- the above framework every spinor is transformed out of the two 
constant basic spinors Ua(C)5 of opposite chirahty A = ±, as follows 

= ^u_,(C), u+(C)=ifu_(C), v' = -i, {vO = o. (1) 

The usual relations hold: (uxiC) = 0, uJxUxiC) = Ua(C), Ua(C)ua(C) =X^a, ^uxip) = 0, 
'^xuxip) = ux{p), ux{p)ux{p) =i>'^\-, where ujx = |(1 + A75). Spinors for the massive 
particle with four- momentum p (with = m^) and spin projection A/2 are defined in 
terms of massless spinors 

TTl Tfl 

where p(^ = p = p — C, m'^/{2(p) is the light-cone projection (p^ = 0) of the p obtained 
with the help of the constant auxiliary vector (. 

The above definition is supplemented in ref. fl^ with the precise prescription on spin 
quantization axes, translation from spin amplitudes to density matrices (also in vector 
notation) and the methodology of connecting production and decay for unstable fermions. 
We collectively call these rules global positioning of spin (GPS). Thanks to these we are 
able to easily introduce polarizations for beams and implement polarization effects for 
final fermion decays (of r leptons, t-quarks), for the first time also in the presence of 
emission of many ISR and FSR photons! 

The GPS rules determining spin quantization frame for u{p, ±) and v{p, ±) of eq. (§) 
are summarized as follows: (a) In the rest frame of the fermion, take the 2;-axis along 
(b) Place the x-axis in the plane defined by the 2;-axis from the previous point and 
the vector ff, in the same half-plane as "ff. (c) With the y-axis, complete the right-handed 
system of coordinates. The rest frame defined in this way we call the GPS frame of the 
particular fermion. See ref. |Ty] for more details. In the following we shall assume that 
polarization vectors of beams and of outgoing fermions are defined in their corresponding 
GPS frames. 

The inner product of the two massless spinors is defined as follows 

S+{PUP2) = U+{pi)u^{p2), S-(P1,P2) = U^{pi)u+{p2) = " (s+ (pi , P2) )* • (3) 

The above inner product can be evaluated using the Kleiss-Stirling expression 

s+(p, q) = 2 (2pC)-'/' {2qO-'/' [(K)(w) - {pv){qO - te,.,X'vW] (4) 

in any reference frame. In particular, in the laboratory frame we typically use ( = 
(1, 1, 0, 0) and r] = (0, 0, 1, 0), which leads to the following "massless" inner product 

s+{p,q) = -{q^ + ^q^)^/{p^ - P^) / {q^ - Q^) + {p^ + ^P^) V(g° - q^) / {p^ - P^) ■ (5) 
Equation (^ immediately provides us also with the inner product for massive spinors 

m(Pi, Ai)m(p2, Aa) = S{pi,mi,Xi,p2,m2,X2), 
m(pi, Ai)t;(p2,A2) = S'(pi,mi,Ai,p2,-"^2,-A2), 
v{pi, Xi)u{p2, X2) = S{pi,-mi,-Xi,p2,m2,X2), 
v{pi, Xi)v{p2, X2) = S{pi, -nil, -Xi,p2, -m2, -A2), 
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where 

In our spinor algebra we shall exploit the completeness relations 

^ + m = u{p, X)u{p, A), ^ — m = f (p, X)v{p, A), 

A A 

f^ = j2u{k,x)u{k,x), e = o. 



(8) 



For a circularly polarized photon with four- momentum k and helicity a = ±1 we adopt 
the KS choice (see also ref. [Q) of polarization vectoiQ 

^'^^^^^ V2u,Ak)nAf3y ^^^^^^ V2u^AkKiCy 

where (3 is an arbitrary light-like four- vector = 0. The second choice with Uo-(C) (not 
exploited in P]) often leads to simplifications in the resulting photon emission amplitudes. 
Using the Chisholm identity^ 

u^{k)-f^u^{P) = 2u^{P) u^{k) + 2u-^{k) u^.{(3), (10) 
M^(A;)7^u<,(C) 1^ = 2u.(C) u,{k) - 2u^^{k) u_.(C), (11) 

we get two useful expressions, equivalent to eq. (0): 



iak, or = MOuM - u^Mu^^io] ■ 



(12) 



In the evaluation of photon emission spin amplitudes we shall use the following im- 
portant building block - the elements of the "transition matrices" U and V defined as 
follows 

w(pi,Ai) ^Ak,P)u{p2,X2) =U = Ul^x,{k,Pi,m,,p2,m2), 

v{p„X,) C{k,0v{p2,X2)=V{'^)[l[Z]=V,l,^{k,p„m„p2,m2). 



^ Contrary to other papers on Weyl spinor techniques |g, ^ we keep here the expHcitly complex 
conjugation in e. This conjugation is cancelled by another conjugation following from Feynman rules, 
but only for outgoing photons, not for beam photon, as in the Compton process, see ref. jl^ . 

^ For /3 = C the identity is slightly different because of the additional minus sign in the "line-reversal" 
rule, i.e. Ua{k)j^^Ua{C) = — u_cr(C)7^u-(T(fc), in contrast to the usual u„{k)^^^u„{f3) — +u^„{(3)^^^u-a{k). 
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In the case of Uo-(C) the above transition matrices are rather simple^: 
U'^{k,pi,mi,p2,m2) = V2 



2Ck 

Ux,,x^{k,pi,mi,p2,m2) = [-U^^^^^{k,p2,m2,pi,mi)Y , 
^a'i,A2(^'P1'"^1'P2,?«2) = t^-Ai,-A2(^>Pi> -mi,p2, -ma). 

The more general case with Ua{(3) looks a little bit more complicated: 

V2 



(14) 

(15) 
(16) 



U^{k,pi,mi,p2,m2) 



X 



s+{pi, k)s_{p,p2) + mim2-^g-g^, miy^s+(A;,p2) + m2y^s+(pi, /c) 



"^i'\/S"*-*^^'^2) +m2 Jgs_(pi,/?), s_(pi,/5)s+(A;,]32) +mim2 



2C/3 2Cfc 
2Cpi 2Cp2 



(17) 



with the same relations (|T^) and (p!6|). In the above the following numbering of elements 
in matrices U and V is adopted 



{(Ai,A2)} 



(++) (+-) 
(-+) (— ) 



;i8) 



When analysing the soft real photon limit we shall exploit the following important diag- 
onality property^ 



u_^{k)u„{Q 

which also holds in the general case of u„{(3), where 

V2 



2Cp 
2Ck 



Sa{k,p), 



(19) 
(20) 



ba{k,p) 



s-a{k,l3) 



s_^{p,p)s,{p,k) + ^^{2pc) m) 



(21) 



3 Born spin amplitudes 

Let us calculate lowest order spin amplitudes for e~(pi)e+(p2) /(P3)/(P4)- For the 
moment we require / 7^ e. Using our basic massive spinors of eq. (^ with definite GPS 

^ Our U and V matrices are not the same as the M-matrices of ref. ^ , but rather products of several 
of those. 

^ Let us also keep in mind the relation b^„{k,p) = —(bcr{k,p))*, which can save time in the numerical 
calculations. 
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helicities and Feynman rules, we define 



^^^^^ X2 - Ms^ + tYnXyMB 



(22) 



where are the usual chiral (A = +1, —1 = R, L) coupling constants of the vector 
boson = 7, Z to fermion / in units of the elementary charge e. 

Spinor products are reorganized with the help of the Chisholm identity (p^ , which 
applies assuming that electron spinors are massless, and the inner product of eq. (|^): 

r r e,B f,B rp rp/ I e,B f,B jji jj 1 

<B PI fx) = 2ie^ \^ ^Ai,-A2[5'ai J-x-jXi ^ 9\i ^AsAz ^AiA4J 

^^^^ ' - Mb" + iTBXyMB ' ^ ' 

where 

^AaAi =U(P3, A3)u(pi, Ai) = 5(^3, mg, A3, Pi, 0, Ai) , 

^A2A4 =v{P2A2)v{Pa,\a) = ^(P2,0, -A2,P4,-m4,-A4), ^^4) 

^A3A2 =w(P3,A3)i;(p2,-A2) = 5'(p3,m3, A3,P2,0, A2), 

f/AiA4 =u{pi,-Xl)v{p4,\4) = 5'(pi,0, -Ai,P4, -"^4, -A4). 

We understand that the total s-channel four-momentum X is always the four-vector that 
enters the s-channel vector boson propagators. Let us stress that the above Born spin 
amplitudes will be used for pi, which do not necessarily obey the four-momentum conser- 
vation P1+P2 = P3 +P4- This is necessary because, in the presence of the bremsstrahlung 
photons, the relation X = pi + p2 = ps + Pa does not usually hold. Furthermore, any of 
the Pi may, and occasionally will, be replaced by the momentum k of one of photons. In 
this case, the spinor into which k enters as an argument is understood to be massless. 

4 First order, one virtual photon 

The 0{a^) contribution with one virtual and zero real photon reads 

[I] (X) = (X) [l + QlF,{s,m,)+Q}F,{s,m,)] + A^box K] W, (25) 

where Fi is the standard electric form-factor regularized with photon mass. We omit, for 
the moment, the magnetic form-factor F2; this is justified for light final fermions. It will 
be restored in the future. In Fi we keep the exact final fermion mass. 

In the present work we use spin amplitudes for 7-7 and 7-Z boxes in the small mass 
approximation m^/s 0,mj/s 0, following refs. [17, 



e,B f,B rp rp/ I e,B f,B jji jj 
KA \PMV\ 0„-^2 5'Ai fi'-Ai ^A3AiJa2A4 + fi'Ai 5'Ai ^A3A2^AiA4 . . 

-Mbox U \X) = lie > — 2 I -r v2ii\/r c)Ai,-A2C'A3 -A4 

fit^z X^-Mb +iVBXyMB (26) 



Ct — — 

-QeQf [S\i,X3 /BDp(M|,m^,S,t,M) - (5ai -A3 /BDp(M|,m^,s,M,t)] 
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where M| = M| — iMzTz, = mi, and the function /bdp is defined in eq. (11) of 



ref. [0. The Mandelstam variables s,t and u are defined as usual. Since in the rest of 
our calculation we do not use rn?j/s — > 0, we therefore intend to replace the above box 
spin amplitudes with the finite-mass results. (NB: For the 7-7 box the spin amplitudes 
with the exact final fermion mass[| were given in ref. ||T^ .) 

5 First order 1-photon, ISR alone 

In order to introduce the notation gradually, let us first consider the 1-photon emission 
matrix element separately for ISR. The first order, 1-photon, ISR matrix element from 
the Feynman rule reads 

viP2, X2) Ml (#1 + m - ^) Cik) u{pu Ai) 
eO ^''^ 

where Mi is the annihilation scattering spinor matrix (including final state spinors). The 
above expression we split into soft IR parts0 proportional to ( ± m) and non-IR parts 
proportional to }i. Employing the completeness relations of eq. ( |l3l) to those parts we 
obtain: 



Mf"" =^ u it) - ^ (^) K?] Kf] 



'2 



P 



(28) 



P P 



where Q5i [^^^^2] ~ "^(Pa, A2)Mim(pi, Ai). The summation in the first two terms gets elim- 
inated due to the diagonality property of U and V , see eq. (pJ]), and leads to 

Gf,^) = KfJ + r(^) gf^a (A:), 



Ikp^ ^ ^ va, LPAij 2fcp2 ^ '-^''^^ '^''^^ (29) 

The soft part is now clearly separated and the remaining non-IR part, necessary for the 
CEEX, is obtained. The case of final state one real photon emission can be analysed in a 
similar way. 



^ It seems, however, that the 7-Z box for the heavy fermion is missing in the Uterature. 
® This kind of separation was already exploited in refs. We thank E. Richter-W§,s for attracting 
our attention to this method. 
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6 First order 1-photon ISR+FSR 

The first order, ISR+FSR, 1-plioton matrix element, with explicit split into IR and non-IR 

parts, reads 

-^S'^ (It) = ^ [I] {P-k)+ ® 11] (P) + r^'^ K] {P-k)+ K] (P), 

(30) 

where we use the compact notation [^] = [^j^^AsaI] ' lowest order Born spin 

amplitudes 03 are defined in eq. (0). The other ingredients are the initial state non-IR 
part: 

p p 

(31) 

and the final state non-IR part 



(32) 



2kp3^ V-/LPA3J LA1A2PA4J V / 2A;p4 



The FSR s-factor 

we define analogously to the ISR case. 

T Coherent exclusive exponentiation, zero and first order 

Spin amplitudes in the zero-th order coherent exclusive exponentiation, C'(tt°)cEBX! 
define as follows 

(34) 

where the s-channel four-momentum in the resonance propagator is = Pi + P2 — 
Yl^=i pi ki- The partition p is defined as a vector (pi, p2, ■ ■ . , pn) where pi = I for ISR 
and pi = for FSR photon, see the analogous construction in refs. For a given 

partition is therefore the total incoming four-momentum minus four-momenta of ISR 
photons. The coherent sum is taken over set {p} of all 2" partitions - this set is explicitly 
the following 

{p} = {(0,0,0,...,0), (1,0,0,. ..,0), (0,1,0,...,0), (1,1,0,...,0), ...(1,1,1,...,1)}. 

(35) 



In eq. (|3^) we profit from the Yennie-Frautschi-Suura m fundamental proof of factoriza- 
tion of all virtual IR corrections in the form-factoi|^ exp(ai?4), where 

54(^1, ...,P4) = QlB2{puP2) + Q}B2{P3,P4) 

+ QeQfB2ipi,P3) + QeQfB2ip2,P4) - QeQfB2{puPi) - Q eQ f B2ip2, Ps) ■ 

9 (36) 

, , r d^k i / 2p + k 2q-k 



A;^ — + (27r)'^ yk"^ + 2kp + ie k"^ — 2kq + ie 

In the above we assume that IR singularities are regularized with a finite photon mass rrij 
which enters into all B2^s and implicitly into s-factors (and the real photon phase space 
integrals, in the following discussion). 

The auxiliary factor F = X^/{p3 + p^)"^ is, from the formal point of view, not really 
necessary. Note that the F-f actor does not affect the soft limit; it really matters if at least 
one very hard FSR photon is present. However, the F-factor is very useful, because it is 
present in the photon emission matrix element, both in 0{a^) and also in all orders in the 
leading logarithmic (LL) approximation. It has also been present for a long time now in 
the "crude distribution" in the YFS-type Monte Carlo generators, see for instance ref . . 
It is therefore natural to include it already in the 0{a^) exponentiation. Otherwise, this 
F-factor will be included order by order. However, in such a case, the convergence of 
perturbative expansion will be deteriorated. As we shall see below, the introduction of 
the F-factor will slightly complicate the first order exponentiation. 

The complete set of spin amplitudes for emission of n photons we define in 0{a^)^^j^y- 
as follows: 

TvfCl) (Pkik2 kn\ _ aB4{pi,...,p4,) 
•^''-n \\aicT2 ■ ■ ■ (Jn) ^ 

n / n 

X E n ^ K] {X,) (1 + 5« ,) + O^Bo. [I] {X,) + \^ (X 

{p} »=i ^ i=i 



s^{k) 



(37) 



The IR-finite 6yl^^ and 3?box are determined unambiguously by identifying for n = the 
above equation with eq. (p5D, up to terms of 0{a^). We obtain 

5y-rt(s) = QlFi{s,my) +(5/Fi(s,m^) - QlaB2{s,m^) - QjaB2{s,mj). (38) 

The 3^Box is obtained from A^box by means of the substitution[^ 

/bdp(M|, m^, s, t, u) /bdp(M|, "^7, s, t, u) - fm{m^, t, u), (39) 



^ In the LL approximation it is, of course, the doubly-logarithmic Sudakov form-factor. 
In the above procedure of subtracting IR divergences, there is no reference to cut on photon energy, 
only reference to the photon mass, similar to the YFS exponentiation on squared spin-summed amplitudes. 



8 



where 



2 2 f t\ f \ 1 f t\ 

fm{m^, t, u) = -B2{m^, t) - -B^^m^, m) = In - In + - In - . (40) 

Similarly the IR- finite Jlf^^ is determined uniquely by identifying, for n = 1, eq. (|37| ) with 
eq. (|30|). In particular the factor F — 1 = (ps +^4 + ujkjY/ {p^ +^4)^ — 1 is a consequence 
of the introduction of the F-factor in eq. (^If). If it was not included, then the 1-photon 
part in eq. ( pTj) would not reduce to the amplitude of eq. (^). Thanks to the presence of 
F — 1, for n = 1, we recover in eq. (pT]) the correct first order amplitude of eq. (p^). 

For very narrow resonances the photon emission in the decay process is separated from 
the photon emission in the production process by very large time-space distance. The 
ISR*FSR interference is therefore strongly suppressed, typically by T/M factors. Since 
our real photons are present down to arbitrarily low k'^^^ = e^fsjl <C F, the effects due 
to the resonance complex phase in the emission of the real photons are taken into account 
numerically and exactly. For virtual photons we have to sum up analytically certain subset 
of the ISR*FSR interferences to infinite order following Greco et al. In practice the 

rule is: multiply each part of the spin amplitude proportional to Z-propagator by the 
additional factor exp{6G{s,t,u)) where: 

6G{s,t,u) = -2Q,Qf-\n i^-j In {^-^ J (41) 

In 0{a^) the above exponential factor induces the additional subtraction in the 7-Z box: 
3Vtbox('5, t, u) Mbox('5, t, u) — 5g{s, t, u). Strictly speaking the above improvement is not 
really necessary, because we would have obtained it order-by-order, through higher order 
virtual non-IR correction. In practice, however, it is mandatory. If we had not made it, 
then the ISR*FSR interference contribution to ApB at Z peak from C(a^)cEEx would be 
dramatically wrong, i.e. 0.5% instead of 0.05%! 

8 Differential cross sections and the YFS form-factor 

The master formula for the unpolarized C'(tt'^)cEEx total cross section is given by the 
standard quantum-mechanical expression of the type "matrix element squared modulus 
times phase space" (contrary to typical "parton shower" approach) 



00 



-If 1 

n=0 A,CTl,...,(Tn = ± 

where the Lorentz invariant phase space (LIPS) integration element is 



„ / n \ n 

dT^{P;p^,P2, ...pn) = j (27r)^5(^) [P - Yj>ij n 



2p0(27r)3- 



(43) 
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The above total cross section is perfectly IR-finite, as can be checked with a little bit of 
effort by analytical partial differentiation]^ with respect the photon mass 



^ (T« = 0. (44) 



dm-y 



Furthermore, the integral of eq. (^21) is perfectly implementable in the Monte Carlo form, 
using a method very similar to those in ref. Traditionally, however, the lower boundary 
on the real soft photons is defined using the energy cut condition > ey/s/2 in the 
laboratory frame. The practical advantage of such a cut is the lower photon multiplicity 
in the MC simulation, and consequently a faster computer program^. If the above energy 
cut on the photon energy is adopted, then the real soft-photon integral between the lower 
LIPS boundary defined by and that defined by e can be evaluated by hand and 
summed up rigorously (the only approximation is m^/nie 0) into an additional overall 
factor exp(2ai?4(pi, ...,^4)), where 

54(^1, ...,Pi) = QlB2{puP2) + Q/^2(P3,P4) 

+ QeQfB2{pi,P3) + QeQfB2{p2,P4) " QeQ fB2{pi, P4) - QeQ/-B2 (P2, Ps) , 



B2ip,q) 



d^k (-1) f p q^ ' 



Stt^ \kp kq 

k0<ey^/2 



Let us introduce 9Jln = e "^"Mn (without virtual IR singularities) and, altogether, the 
above reorganization yields the new expression for the unpolarized total cross-section 



00 

— n ' ^ \ L 



71=0 X.cr 



(46) 



where Y{pi, ...,^4) = 2a_B4(pi, ...,^4) + 2a3fJi?4(pi, ...,^4) is the conventional YFS form- 
factor defined analytically in terms of logs and Spence functions - we do not show it here 
explicitly due to lack of space, see refs. [|^|1T],^, In the YFS form-factors we keep 
the final fermion mass exact. The fully exclusive differential cross section of eq. (|46|) is 
already implemented in the Monte Carlo event generator /C/C [|rT| . 

The extension of the above exponentiation procedure to 0{a'^)^j^j^^ and beyond re- 
quires more work, but does not pose any conceptual problem. It will be implemented in 
the future version of the JCJC Monte Carlo. 



This method of vahdating IR-finiteness was noticed by G. Burgers The classical method of 

ref. relies on the techniques of the Melin transform. 

^'^ The disadvantage of the cut fc" > ey/s/2 is that in the MC it has to be implemented in different 
reference frames for ISR and for FSR - this costs the additional delicate procedure of bringing these two 
boundaries together, see ref. |]ll| and/or discussion in the analogous i-channel case in ref. 



10 



9 Fermion spin polarization and photon spin randomization 



The great advantage of working with spin amphtudes is the easiness of introduction of 
full spin polarizations for all particles. The general case of the total cross section with 
polarized beams and decays of unstable final fermion being sensitive to spin polarization 
Pl|23|-E5l reads 



n=0 (Ti a,b,c,d=0 X^^\- 

^a^b b _ m>(r-) (pkik2 A;„\ raT>(r) (Pkik2 icid 



(47) 



where, for k = 1,2,3, a'' are Pauli matrices and 0"°^ = 6x^^ is the unit matrix. The 
components e", a, & = 1,2,3 are the components of the conventional spin polarization 
vectors of e~ and respectively, defined in the so-called GPS fermion rest frames (see 
ref [|T0| for the exact definition of these frames). We define e° = 1 in a non-standard way 
(i.e. Pi ■ ii = rrie). The polarimeter vectors hi are similarly defined in the appropriate 
GPS rest frames of the final unstable fermions {pi - hi = rrif). Note that, in general, hi 
may depend in a non-trivial way on momenta of all decay products, see refs. p^,^ for 
details. We did not introduce polarimeter vectors for bremsstrahlung photons, i.e. we 
take advantage of the fact that the high energy experiment is completely blind to photon 
spin polarizations. 

Let us finally touch briefiy upon one very serious problem and its solution. In eq. ( ^T] ) 
the single spin amplitude Wln^ already contains 2"'{n + 1) terms (due to 2" ISR-FSR 
partitions). The grand sum over spins in eq. ( ^Tf ) counts 2"4'*4^ = 2""'"^^ terms! Altogether 
we expect up to ~ n2'^^~^^^ operations in the CPU time expensive complex (16 bytes) 
arithmetics. Typically in e~e^ — fi'^ the average photon multiplicity with > IMeV 
is about 3, corresponding to ~ 10'^ terms. In a sample of 10^ MC events there will be a 
couple of events with n = 10 and A^ = 10^^ terms, clearly something that would "choke" 
completely any modern, fast workstation. There are several simple tricks that help to 
soften the problem; for instance, objects such as J2a ^i'^xx 3-factors are evaluated 

only once and stored for multiple use. This is however not sufficient. What really helps 
to substantially speed up the numerical calculation in the Monte Carlo program is the 
following trick of photon spin randomization. Instead of evaluating the sum over photon 
spins CTj, i = 1, ...,n in eq. (^Tj), we generate randomly one spin sequence of (ai, ...,an) 
per MC event and the MC weight is calculated only for this particular spin sequence! 
In this way we save one hefty 2" factor in the calculation timeQ. Mathematically this 
method is correct, i.e. the resulting cross section and all MC distribution will be the 
same as if we had used in the MC weight the original eq. ( ^Tj ) (see a formal proof of the 
above statement in Sect. 4 of ref. ||2^). Let us stress again that it is possible to apply 
this photon spin randomization trick because (a) the typical high energy experiment is 

^•^The other 2" factor due to coherent summation over partitions cannot be ehminated, unless we give 
up on narrow resonances. 
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blind to photon spin polarization, so that we did not need to introduce in eq. ( ^7|) the 
polarimeter vectors for photons, and (b) For our choice of photon spin polarizations the 
cross section is rather weakly sensitive to them, so the method does not lead to significant 
loss in the MC efficiency. 

10 Conclusions 

We presented the first order coherent exclusive exponentiation CEEX scheme, with the full 
control over spin polarization for all fermions. This new method of exponentiation is very 
general and has many immediate and longer term advantages. The immediate profit will 
be the inclusion of the ISR-FSR interferences and availability of the exact distributions 
for multiple hard photons without giving up on exclusive, YFS-style, exponentiation. In 
particular it is applicable to diffucult case of the narrow resonances. The resulting spin 
amplitudes and the differential distributions are readily implemented in the MC event 
generator. (Numerical results will be presented elsewhere.) 
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